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Abstract 



o 
o . 

f^ , We study the effect of magnetic field on the bulk viscosity of nuclear 

ri ^ matter in neutron stars. We employ the framework of relativistic mean 



field theory to observe the dense nuclear matter in neutron stars. The 
effects are first studied for the case when the magnetic field does not exceed 



o 

C/3 . the critical value to confine the electrons to the lowest Landau levels. We 

then consider the case of intense magnetic field to evaluate viscosity for 
• '^ ■ the URCA processes and show that the inequality PF{e) +Pf(p) > PF{n) 

t , ■ is no longer required to be satisfied for the URCA processes to proceed. 



1 Introduction 

After the discovery of niagnetars there has been a lot of interest in the study 
of the effects of strong magnetic fields on the various astrophysical phenom- 
ena. The observed soft gamma repeaters discovered recently by Kuoveliotou C 
et al [1] and Kulkarni and Thomson [2] show the presence of strong magnetic 
fields of upto 10^^ Gauss. Though the core of neutron stars is poorly under- 
stood, magnetic fields as strong as i? ^ 10^^ Gauss may be assumed in the 
core by the scalar virial theorem. It has been argued that even higher fields ( 
B ^ 10^" Gauss) could be consistent with the virial theorem [3]. It is there- 
fore important to investigate the global properties of these objects, particularly 
a detailed investigation of the stability of dense stellar matter in the presence 
of ultra strong magnetic fields. This requires a study of the bulk viscosity of 
such objects, since the damping of vibrations and secular instability of rapidly 
rotating neutron stars depend upon the bulk viscosity. The bulk viscosity of 
neutron stars arises mainly due to deviation from /? equilibrium and the ensu- 
ing non-equilibrium reactions. These important reactions are so called URCA 
processes: 



n — > p + e + Ue (1) 

p + e^ — > n + i^e (2) 

and the modified URCA processes: 

n + n ^ n + p + e~ + i7e (3) 

n + p + e^ -^ n + n + i/e (4) 

At low temperatures for degenerate nuclear matter, the direct URCA process 



can take place only near the fernii energies of participating particles. Simulta- 
neous conservation of energy and momentum puts the constraint 

PF{e) + Pf{p) > PF (n) (5) 

which must be satisfied for the URCA processes to take place. In most 
models of nuclear matter this constraint is satisfied only at high densities, viz, 
Ub > 1.5 no, where uq ~ 0.16 /?n^'^,is the nuclear saturation density. There 
have been many studies of the bulk viscosity of neutron stars under the usual 
situation of zero magnetic field [4-8]. In this note we report the calculations 
of bulk viscosity of neutron stars from the URCA processes in the presence 
of magnetic field. Such a study is of interest because significant changes in 
the equation of state of a magnetized star could lead to similar changes in the 
bulk viscosity as well. We have considered a wide range of magnetic fields in 
our study, from "low" to "very high" magnetic fields (> lO^^G). As has been 
discussed above, such extremely strong magnetic fields in the interior of stellar 
objects can not be completely discounted and have infact been already studied. 
For the description of nuclear matter at high densities which are found in the 
interior of neutron stars, we use the frame work of relativistic nuclear mean field 
theory, including the p contribution, as extended to include strong interactions 
among the particle species involved. As is well known, this would give a fairly 
good description of strong interactions among the various particles and can be 
easily extended to include higher resonances as well. In Section 2 we derive 
expressions for bulk viscosity both in the high and low magnetic field limits. In 
Section 3 we present our results followed by a brief discussions. 

2 Derivation of Bulk Viscosity 

In this Section we derive expressions for the bulk viscosity in the presence of 
magnetic field. The presence of magnetic field affects the bulk viscosity only 
through the changes in the reaction rates as well as through the modification 



of the chemical composition of neutron star matter. We therefore adopt an 
approach similar to that of Wang and Lu B, Sawyer J5|, MadsenS and Gupta 
et al H] for the zero magnetic field case. The bulk viscosity ( is defined by 



J, )av \ A ) 2 

at vq Av oj^ 



C = 2(^U-(^f-, (6) 



Here vq is the specific volume of the star in equilibrium configuration, Av is 
the amplitude of the periodic perturbation with period r = — : 

2'Kt 

v{t) ^vo + Av sin( ) (7) 

r 

The quantity {-^)av is the mean dissipation rate of energy per unit mass 
and is given by the equation 

,dW . 1 f ^, sdv , , , 

i-^u---Jm^,dt (8) 

In the above the pressure P(i) can be expressed near its equilibrium value 
Po, as 

dP dP dP dP 

Pit) = Po + i^)oSv + i — )oSn, + i — )oSn. + (— )o^n„ (9) 

The change in the number of neutrons, protons and electrons per unit mass 
over a time interval (0,t) due to the URCA reactions (1) and (2) is given by 

— Sun — Srip = Sn^ = / —r~dt (10) 

Jo dt 

The net rate of production of protons, —jf- is given by the difference between 
the rates P^ and P^ of reactions (1) and (2) respectively. At equilibrium the 
two rates are obviously equal and the chemical potentials satisfy the equality 

A^ == /i„ - /Ip - /ie = (11) 

A small volume perturbation brings about a small change in the chemical 
potentials and the above inequality is no longer satisfied, now Afi is not zero 
and consequently the reaction rates of the processes (1) and (2) are no longer 
equal. The net rate of production of protons will thus depend upon the value 



of Afi. In the linear approximation, -^^ <C 1, the net rate is proportional to A/i 



feT 

and thus can be written as 



^ = A - A = -\A^, (12) 



(13) 



Using the thermodynamic relation 

dP _ o>M» 
drii dv 

in equation (9) and employing equation (10), we obtain 
SP = P{t)-Po 

^-^J\s,m (14) 

The change in the chemical potential Sfi{t) arises due to a change in the 
specific volume Sv and changes in the concentrations of various species, viz, 
neutrons, protons and electrons. Thus 

Mi) = 5f^{0) + AoSv + Ao6n^ + AoSn, + AoSnAl5) 

OV OUn OHp one 

Using equations (10)-(12), we arrive at the following equations for S^: 

^ ^ ujA— cosiujt) - C\S^ (16) 

at vq 

where 

A-M'-^k (17) 

C^v(— + — -^^) (18) 

dUp dUe dUn 

Since for small perturbations. A, A and C are constants, equatioin (15) can 
be solved analytically to give 

Sfx ^ „ '^t,.9 — [-CAe-^^* + u; sin(cjt) + CX cos{u;t)] (19) 

OJ'' + C'^A^ Vq 

Finally using equations (6), (8), (10), (11), (14) and (17)-(19), we obtain the 
following expressions for C 

^^A LuCX 1 - e-^^^ 



Given the number densities of the these particle species in terms of their 
respective chemical potentials, one can determine the coefficients A and C; given 
the rates Fi and r2 for the two processes (1) and (2) one can determine A and 
hence C for any given baryon density and temperature. For the description of 
the dense nuclear matter we have used Walecka's mean field theory [8-10]. In 
the mean field approximation we allow the mesons fields a, uj and p to acquire 
density dependent average values. From the effective langrangian one can then 
read off the effective masses and chemical potentials ffii and fli. The fermi 
momenta, are in turn related to these effective quantities by 



/2f-mf (21) 

The various chemical potentials and masses required in these expressions are 
obtained from the self consistent solution of the mean field equations: 

TO^CT + 5m„6^^CT^ + cg^fqd-^ ~ gaN{nl + n''^) = 

^Ipo = -^9pN{np - n„) 

rrin - fhn ^ mp - fhp ^ gcm^o 

1 

A^n — A*" ^ Qi^N^a + t:9pNPo 

p.p = fip- gcjNi^a - -^gpNpo 

(22) 

The number densities rii and the scalar densities nf |9[| are given in terms of 
these quantities by the expressions 



-.-2/^(l + e^)- (23) 



("2^ 



"^^-2'^'/(S^i(^+^'^)"' (^^) 



where 



E, = . /p2 + ^2 (25) 



Since the chemical potentials fi„ and fip are of the order of a few hundred 
MeV, even for temperatures upto a few MeV the matter is completely degener- 
ate, this leads to the following expressions for rii and n| : 

1 3 

"'==^(^»-™^)' (26) 

n^ = ^Mf^~rr^)^ (27) 

where i=n, p. To evaluate rie, set fie — Me ^nd rFig — rrie in equation(26). 

It now remains to introduce the magnetic field. In the presence of constant 
magnetic field B in the z-direction, the energy of a charged particle of mass m 
and charge e is given by 

E^ ^pl+m^ + 2veB (28) 

where the quantum number v is given by 

v = VL + ^ + (y (29) 

for the Landau level v^ — 0, 1, 2.... and spin a = ±i. 

In the presence of a magnetic field, the number densities of electrons and 
protons are given by 

"- = l5 E (2 - S,^o)Vl^l -ml- 2veB (30) 

and 

"P = y4 E (2 - 5,.fl)\jnl -rhl- 2ueB (31) 

whereas the neutron number density expression remains unchanged namely 

1 
37r 

In the decay rates the usual sum over the charged particle states (per unit 

volume) in zero field ,viz -p- / <Pp is replaced by Y^ui'^^^^fi) J^^^ dPz J JhL^ dpy 

2 

and the matrix elements have to be calculated by using the exact wave functions 
for electrons and protons obtained by solving the Dirac equation in the presence 
of magnetic field [11-15] . We are now in a position to evaluate the various 
qualities required for the calculation of bulk viscosity, viz. A, C and A. 



^iK-Kr^ (32) 
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2.1 Weak Magnetic Field 

For weak magnetic field several Landau levels are populated and the matrix 
elements remain essentially unchanged [13-14] and one needs to account for the 
correct phase space factor. For non-relativistic degenerate nuclcons the decay 
rate constant A is given by 



17 

EP-M-.^ \ , „ (33) 

which reduces to well known expressions [16-17] in the limit i? ^ 0. Using 
equations (30), (31) and (32) the coefficients A and C given in equations (17) 
and (18) become 



^ _ Jill-fhl) , lZZl^i'^-5ufi)^lll-ml-2ueB 



E::=r(2-5.,o)Me(M?-m2-2z.eB)-^^ 2^^^^'' 



C ^ 



2^2 



M^Mn-^n ^^ E:^=r(2 - &.fi)MPy -H- ^'^eByi 



+ (2£fl)^ (35) 

m,- 



"p 



eB y:ZI^{2 - <5.,o)Me(M? -ml- 2veB)-"^ 
2.2 Strong Magnetic Field 

For strong magnetic field , the electron is forced into the lowest Landau level 
corresponding to i^ = . This happens for a field B > Be , where Be is given by 

2_ 2 

eBc = ^° „™° . Because of the equality np — Ue , dictated by charge neutrality, it 
is then straight forwad to verify that for B > Be only the lowest Landau level u = 
contributes for the protons as well. This leads to considerable simplification 
in the expression for the decay rate. Using the exact wave functions for protons 



and electrons in the lowest Landau level and following [15] , the rate for the 
process (1) can be calculated and we get 

r = ^^-^eB^yf / dEpdEedEndE,ElS{Ep + E, - E„ - E,){1 - fn)fpfe 
{2Tr)° PfKP) J 

[4.9i + {9A + 5.)'][e-''"'"-«'"''%(p|'N - 4p|,(e)) + e-^0(pf (n))] (36) 

Carrying out the energy integrals for degenerate matter , the decay constant A 
is given by 

(p|,(,i)-4p2 (e)) p2(„) 

[e ^ e{pl{n) - Aplie)) + e'^^ e{pp{n))] (37) 

The constants A and C are obtained from equations (34) and (35) by setting 
V ~ Q. Substituting A, C and A into the expression (20) the bulk viscosity can 
now be calculated for a given baryon density and temperature. 



3 Results and Discussions 

In this paper we have used the relativistic mean field approximation to describe 
nuclear matter at high densities. This approach gives a good description of 
nuclear matter and neutron star properties at high densities. We have used the 
values of the various couplings g^^ , QuN , 9pN , b and c given in [10] namely 

^^ = 0.0152502Mey"i 
9ujN 



O.OUMeV-^ 

i£^ = O.OUMeV^^ 
nip 

6 = 3.418 X 10^3 
c = 0.0146 

(38) 

With this set of parameters, one obtains ns = 0.153 fm~^ ; binding energy 
( — — TO„) = —16.3 MeV ; charge symmetry energy =32.5 MeV ; bulk modulous 
=300 MeV and Landau mass = 0.83 to^v. 



We have computed the bulk viscosity for zero magnetic field, for low mag- 
netic fields by summing over various Landau levels and finally for very high 
magnetic fields, where only the ly — level contributes. In figure 1 we show 
the variation of the coefficient of bulk viscosity with baryon density for eB = 
and eB = lO'^ MeV'^ for T = I MeV. As can be seen from (33) and (36) the 
variation of ( with T goes as T^; so for other temperatures C can be obtained 
by scaling. We, in our study, have divided the magnetic field into two regimes, 
the "low" magnetic field where a large number of Landau levels contribute and 
the "strong" magnetic field where only the lowest Landau level contributes. 
The critical field Be separates the two regimes depends crucially on the chem- 
ical potential of the electron/proton. These in turn depend upon the number 
density ns- As ns varies from 0.2 fm^^ to 2.0 fm^^, the range of ub over 
which figure 1 is plotted, eB^ remains well above 10'' MeV'^. Therefore we have 
used the weak field formula for computing viscosity at this value of magnetic 
field. At zero magnetic field the viscosity changes smoothly with ns- As ns 
changes from 0.2 fm^^ to 2.0 fm^^ viscosity changes by nearly one order of 
magnitude. At eB = 10* MeV'^, superimposed on a nearly smooth change is a 
rapid fluctuation in viscosity with hb ■ As one can see, whenever a new Landau 
level appears, the reaction rate shows a sudden discontinous change which is 
reflected in a rapid fluctuation in viscosity. In figure 2, we show the variation 
of viscosity with ng for "high" magnetic fields. This is achieved by considering 
only very small number densities, viz, ns < 0.01 fm^^. For such small hb, 
even at magnetic fields of a few hundreds MeV'^ , only the lowest Landau level 
contributes and therefore one can use the formulas obtained in the high mag- 
netic field case. The direct URCA process is now no longer inhibited whereas 
for B = the process is not allowed by energy momentum conservation and 
the bulk viscosity does not get any contribution. We find that the viscosity 
decreases very rapidly with the changes in number density. For example, for 
eB = 100 MeV'^, the viscosity changes by 12 orders of magnitude as hb varies 
from 0.001 fm^^ to 0.004 fm^^. The variation in viscosity with hb is most 
drastic at lower magnetic fields; at higer magnetic fields the variation with ub 
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becomes much less, so much so that for eB = 1000 MeV'^ the bulk viscosity 
remains almost constant and beyond that, it infact starts increasing with ns- 
In figure 3 we plot bulk viscosity against the magnetic field for two baryon 
densities, viz, ns = 0-34 fm^'^ and ns ~ 0.6 j"mr^ . For magnetic fields upto 
10,000 Me^^, the weak field approximation is valid as pointed out earlier. Here 
we find, not surprisingly, that the viscosity is nearly independent of magnetic 
field. As the magnetic field increases (beyond 10'* MeV^ for ns = 0.34 fm^^ 
and 3 x 10"* MeV"^ for hb — 0.6 fm~^) bulk viscosity begins to increase with 
magnetic field because now only the lowest Landau level contributes and the 
magnetic field can be considered "high" . In figure 4, we show the variation of 
viscosity with magnetic field at low densities. We find that initially the vis- 
cosity increases very rapidly but then saturates to a nearly constant value as 
the magnetic field increases. To conclude we find that in the presence of strong 
magnetic field, the triangular inequality is no longer required to be satisfied for 
the URCA process to contribute to the viscosity. At very low densities this re- 
quirement can be met for magnetic field strengths of a few hundred MeV"^ , and 
the bulk viscosity of the medium is finite though small. This is to be contrasted 
with the situation in which the magnetic field is zero and the URCA process 
does not contribute to bulk viscosity upto a certain minimum density. However 
at densities of relevance in neutron stars (n^ '^ 2 — 4 ng), the effect of non-zero 
magnetic field will be felt only at extremely high values of B(> 10^ ~ 10^). 
even though the reaction rate is directly proportional to the magnetic field, the 
viscosity depends in a rather complicated fashion on the magnetic field. 
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Figure 1 

Figure 1: Bulk Viscosity C in units oi gm.cm.s"^ as a function of Baryon Density 
ub in fm-^ at T = 1 MeV for magnetic field eS = and eB = lO'^ MeF^. 
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Figure 2: Bulk Viscosity C as a function of Baryon Density ng for different 
values of the magnetic field at T = 1 MeV. The units of (, ng and B are same 
as in Fig. f . 
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Figure 3: Bulk Viscosity C as a function of magnetic field at T = 1 MeV for 
different values of baryon density. The units of C, ns and B are same as in Fig. 
1. 
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Figure 4: Bulk Viscosity C as a function of magnetic field at T = 1 MeV for 

different values of baryon density. The units of C, ns and B are same as in Fig. 

1. 

16 



